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ABSTRACT 

The effects on stability of inertial forces arising 
from closed-loop activation of mass-unbalanced con- 
trol surfaces are studied analytically using inertial en- 
ergy approach, similar to the aerodynamic energy ap- 
proach used for flutter suppression. The limitations 
of a single control surface like a leading-edge (L.E.) 
control or a trailing-edge (T.E.) control arc demon- 
strated and compared to the more powerful combined 
L.E.-T.E. mass unbalanced system. It is shown that 
a spanwise section for sensor location can be deter- 
mined which ensures minimum sensitivity to the mode 
shapes of the aircraft. It is shown that an L.E. con- 
trol exhibits compatibility between inertial stabiliza- 
tion and aerodynamic stabilization, and that a T.E. con- 
trol lacks such compatibility. The results of the present 
work should prove valuable, when flutter is suppressed 
using mass unbalanced control surfaces, or for the sta- 
bilization of structural modes of large space structures 
by means of inertial forces. 

NOMENCLATURE 


a a 

scalar gains, see eq. (36) 

a(w) 

transfer function defined in eq. (45) 

b 

reference semichord length 

h 

element ij of the [ B c ] matrix 

t>R 

semichord length at the reference 


section 

by 

n 

c ij 

element i, j of the real part of [T] 

dij 

element i,j of the imaginary part 


of [T] 

11 

scalar gains, see eqs. (36) 

h 

plunge coordinate, positive down 

h-R 

plunge coordinate at the reference 


section, positive down 

~b.R 

= blR/bR 

ky 

- dh 

dy 

M 

mass of the strip 


mass of the L.E. control 


Mr mass of the T.E. control 

Ml = Ml/M 

Mt — Mt/M 

n number of structural modes 

p t transfer function pole location 

r number of control 

s reference length, of the Laplace 

variable 

tti &\\ + iea 

t x j defined in eq. (36) 

y spanwise coordinate, measured from 

the reference section 

yoL spanwise distance of L.E. control 

surface center of gravity 
measured from reference section 

ycj spanwise distance of T.E. control 

surface center of gravity 
measured from reference section 

j/GL = Ugl/K 

VGT = VGT/br 

a pitching angle, positive nose up 



aR angle of attack at reference section 

of activated strip 

Oty = <Xy/otR 

(3 leading-edge control deflection, 

positive nose down. 

8 trailing-edge control deflection, 

positive trailing edge down 

\ eigenvalue 

\ i largest eigenvalue of a binary system 

A 2 smallest eigenvalue of a binary 

system 

p fluid density 

u ) frequency 

u; n natural frequency 



Matrix Notation 

[Ar ], [ Ai ] real and imaginary parts of aero- 
dynamic matrix, respectively, 
(order nx (n+ r)) 

[Ar,s], [Ai_s] real and imaginary parts, respec- 
tively, of the aerodynamic 
matrix associated with the 
structural modes (order nx n) 

[Ar c], [Ai c ] real and imaginary parts, respec- 
tively, of the aerodynamic matrix 
that couples the control surfaces 
with the structural modes 
(order nx r) 

[B] mass matrix (order n x (n+ r)) 

[BJ structural mass matrix (order nx n) 

[B c ] coupling mass matrix, between 

structural modes and control 
surfaces (order nx r) 

[E] stiffness matrix (order nx ( n + r) ) 

[EJ structural stiffness matrix 

(order nx n) 

[E c ] coupling stiffness matrix, between 

structural modes and control 
surfaces (order nx r) 

{F} column matrix of forces (order n) 

[Qfi]. [Q/l real and imaginary parts, respec- 

tively, of the energy eigenvectors 
(of matrix [Q] (order n x n)) 

{g} column of structural responses 

(order n) 

{g} column of structural and control 

responses of order (n + r), 
see eq. (3) 

(go } column of structural response 

amplitudes (order n) 

{g c } column of control surface responses 

(order r) 

[T] matrix representing control law 

transfer functions (order r x n) 

[U] energy eigenvalue matrix, 

see eq. (8) (order n x n) 


{6? },{£/} real and imaginary parts, respec- 

tively, of the energy principal 
coordinates, see eq. (12) 

Superscripts 

T transposed matrix 

\ 

* conjugate 

Subscripts 

c relates to control surfaces 

s relates to the main structure of 

the system 

R relates to the reference section 

within the 2-D strip 

INTRODUCTION 

Studies concerning flutter suppression using ac- 
tive controls concentrate on the numerical design of 
control laws and their effects on the resulting perfor- 
mance of the system (Mukhopadhyay and others 1981, 
Newsom 1979, Newsom and others 1980, Nissim and 
Abel 1978, Mahesh and others 1981, and Freymann 
1987). In these studies, the system is considered to be 
fixed, and no attempt is made to investigate the possi- 
ble effects of varying some of the aircraft parameters 
on the resulting control law and the overall system per- 
formance. Even for existing systems, there are some 
parameters that can bc readily changed, like the mass 
balancing of control surfaces or the control surface ac- 
tuator dynamics. 

In the following, the effects of control surface 
mass unbalance of active control systems on the 
closed-loop stability of the systems are studied. This 
problem is relevant because aerodynamic control sur- 
faces are naturally unbalanced. The balancing of these 
control surfaces requires the addition of non-negligible 
masses. Because an investment is made in an active 
control surface in terms of the mass of the elements 
required to drive the control surface, one is naturally 
tempted to avoid adding further masses to balance the 
active control surface. An analytic approach to the 
problem was adopted so that the results obtained would 
be informative regarding general systems rather than 
for a specific specialized system. For this reason, a 
two-dimensional model is considered, representing a 
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possible strip of a finite wing. To avoid assigning fre- 
quencies or mass values to this strip (since continuous 
wings have, effectively, an infinite number of frequen- 
cies, masses, and mode shapes), an energy approach 
was adopted, whereby the work done by the inertial 
mass-unbalanced forces per one cycle of oscillation 
was evaluated (Nissim 1971 and 1977). If these forces 
do a positive work on their surroundings, the system 
is dissipative. If, on the other hand, these forces do a 
negative work on their surroundings — implying that 
a positive work is done on the system — the system 
absorbs energy and is unstable. Hence, the sign of 
the work done by the unbalanced inertial forces indi- 
cates their potential contribution to either stability or 
to instability. 

The major part of the work presented herein re- 
lates to a 2-D rigid strip performing pure pitch-plunge 
oscillations. In the final part of this work, the strip is 
allowed to twist along its span, to vary its chord length 
along the span, and to rotate as a result of the wing’s 
bending oscillations. It is shown that this more realis- 
tic, generalized strip geometry and motion can be re- 
duced for simple geometries and appropriate selection 
of reference section to yield results identical to the 2-D 
pitch-plunge rigid strip. 

INERTIAL ENERGY APPROACH 

Basic Equations 

Let the n equations 

{F} = -w 2 [ B + upb 4 s( A r + iAi) ] {q} + [ E] {<?} 

( 1 ) 

represent the generalized equations of motion of n 
structural modes with r-activated controls, where 
at flutter 

{F} = 0 

and 


[E] stiffness matrix 

p fluid density 

s reference length 

b reference semichord length 

{?} generalized coordinates 

All matrices in eq. (1) are of order n x ( n+ r) , that is, 
n structural modes and r active controls. Equation (1) 
can therefore be written as 

{F} = (^-u} 2 ^[B s \B c ] + -npb 4 s{[A R<s \A RyC ] 

+ i[Aj tS \Ap c \) 

where 

ts> = | 

Assume a control law of the form 

{< lc]={T}{q } ( 4 ) 

where [T] is a r x n matrix representing the transfer 
functions of die control law. The matrix [T] is there- 
fore, in general, a function of the Laplace variable s. 
By substituting eq. (4) into eq. (2), the following equa- 
tion is obtained: 

{F}= + [B c ][T] + i Tpb 4 s({A R}S } 

+ [A.R'c\[T] + i [Aj tS ] 

+ i[Ar, c )[T}) 

+ IF,] + lF c ][T’]^{g} 

(5) 


+ [ E S \E C ] ) ■] q 

Qc 


( 2 ) 


Q 

Qc 


(3) 


u; 

[B] 

[A/i], [A/] 


frequency 


mass matrix 


real and imaginary parts of aero- 
dynamic matrix, respectively 


It is shown in Nissim (1971 and 1977) that the 
work P done per cycle of oscillation by the system on 
its surroundings is given by 
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P = 5^ U ;j 


— ( [Aj fS ] + [Ai fS ] T + [A IiC ][T] 


+ i 


+ [r] T [A/, c ] 

/ 

/ [Bcim ~[r*] T [g c ] T 

\ 7rp6 4 s 

+ [ A R:S \ — [j4/i iS ] T + [A Rfi ][T} 


[T*] t [A Ri c] T 


{go} 


( 6 ) 


where 

{g} = {goK"‘ (7) 

It is also shown in Nissim (1977) that if [U] is defined 
to represent the expression within the square brackets 
in eq. (6), that is 


[U] 


[Ai tS ] + [ Ai iS ] r + [Ai, c ][T] 


where X,- is the ith eigenvalue of eq. (10), neces- 
sarily real since [U] is Hermitian, and £ R and £i are 
defined by 


{go} = [Q R +iQi}{tR+ *6} ( 12 > 

The columns of the square matrices [QkI and [Q/ ] 
are the real and imaginary parts of the eigenvectors of 
eq. (10). The vector {£ R + i£i} contains the general- 
ized modal coordinates of the transformation defined 
by eq. (12). 

It should be noted that if all the As in eq. (11) are 
positive, the system is stable (dissipative), irrespective 
of the responses represented by the different £ values. 
If, however, one or more of the As assumes a nega- 
tive value, the system may become unstable if the £ 
responses are such that P becomes negative. If one 
wishes to ensure stability irrespective of the responses 
of the system, all the X values must be made positive. 

In the following, it is desired to determine the 
contribution of the out-of-balance inertial terms inde- 
pendent of the aerodynamic terms (which are depen- 
dent on the flight configurations), and, therefore, con- 
sider the following expression for [U]: 


+ i 


+ [T*] r [A IiC ] 

/ 

[ B C ][T ] -[T*] J [B C ] J 
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pb A s 


+ [A r , s ] 


[U] = i[B c ][T]-i{T*] T [B c } J (13) 

where the work done by the inertial out-of-balance 
terms is given by 


— [Ar ^] 1 + \.A R 'c\[T\ 

- [r] T Ufl, c ] T ) 


( 8 ) 


so that 


P = ^P b * u l± Lgo^tt/Hgo} (9) 

then, by taking the eigenvalues of [U], that is, 
by solving 

X{7]} =[(/]{ Tj} (10) 

P can be reduced to the form 

P= ^£[, l(dl + ^) + x 2 (& + &) 

+ . . . + A n (£^ n + £/ n )j (11) 


P=^[Qo\[U]{qo} ( 14 ) 

Object of Analytical Work 

The object of the following work is three-fold: 

1 . To determine a control matrix [T] that will ensure 
stability, that is, so that all the X eigenvalues are 
positive. 

2. To determine the effect of mass unbalance on the 
Xs using the aerodynamic energy transfer func- 
tion matrix [T], 

3. To determine the possible compatibility between 
(1) and (2) above. 
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DETERMINATION OF [T] THAT 
ENSURES STABILITY OF MASS- 
UNBALANCED CONTROL SYSTEM 

Basic Equations 


and also 

C2l = d,2 1 = c 22 = <222 = 0 

Similarly, for the case of a T.E. control only (that is 
when (3=0) 

&ii = 621 = 0 


Assume a two-dimensional strip perlorms pitch- 
plunge oscillations with activated L.E.-T.E. con- 
trol surfaces. 

Let 



where (3 is the L.E. rotation, positive nose down, and 6 
is the T.E. rotation, positive T.E. down. The deforma- 
tions h and a indicate plunge (positive down) and pitch 
(positive nose up, Fig. 1). The parameter b denotes the 
scmichord length of the 2-D wing. 

Denote 


and also 

ci 1 = d\\ = C12 = d \2 = 0 

The characteristic equation for the determination of the 
eigenvalues A of matrix [U], as given by eq. ( 13 ), is 
obtained from 

\\I-U\ = 0 

That is, 

\ 2 + 2 A( 6n d\\ + 612 d 2 i + b2\d\2 + 622^22) 

+ 4 ( 6n dn 621^12 + bll dll *>22^22 

+ bnd2ib2\d\2 + bnd2\b2id22) 

2 

_ (&ndi2 + bi2d22 + 621 dll + 622 d2i) 

— ( 621 CH + & 22 C 21 — b\\C\2 b\2C22) ~ ® 

( 19 ) 


where some detailed derivations for the different 6 , ; 
terms can be found in appendixes A and B. Let [T], for 
any value of frequency, be given by 


m 


ci 1 + idn C 12 + idn 

C21 + id21 C22 + *d22 


(17) 


Substituting eqs. (16) and (17) into cq. (13) one ob- 
tains, after some rearrangements, the following equa- 
tion for [U]: 


2 ( 61 1 dn + &i 2 d 2 i) 


\U] = - I - 1(621 Cll + 622 C 21 - 611 C 12 - bi 2 C 22 ) 
+ (621 dn + 622 d 2 i + b\\dn 
+ 612 d 22 ) 

t( 621 Cll + f>2 2C21 — &11C12 — &I 2 C 22 ) 

+ ( 61 1 di 2 + b 1 2 d 22 + &21 di 1 
+ 622 d 2 i) 


2 ( &21 dl2 + ^22 d22 ) 


(18) 


It should be noted (see eqs. (2) and (15) to (17)) that for 
the case of an L.E. control only (that is when 6 — 0) 


612 = 622 = 0 


which can be rearranged to yield 

A 2 + 2 A( 6n dn + bi 2 d 2 i + b 2 \di 2 + 622 d 2 2) 

+ 4( i>n dn 622 d 2 2 + 6i 2 d2i62idi2) 

— ( 611 di2 - £>21 dll) 2 - (f>12d22 - £>22d 2 l) 2 

- 2(6udi2 + 62idii)(6i 2 d 2 2 + b 2 2 d 2 \) 

- (621 cii + i>22C2i - buen - &12C22) = 0 

(20) 

Equation ( 20 ) forms the basis for the stability investi- 
gations presented in the following sections. 

The Case of T.E. Control Only 

To simplify the study of eq. ( 20 ), consider first the 
case of an activated T.E. only. In this case, 

6n = 621 = cii = C12 = dn = di2 = 0 (21) 

Substituting cq. ( 21 ) into the characteristic eq. ( 20 ), 
one obtains 

A 2 + 2 A( 612 d 2 i + ^22 d 2 2 ) 

— ( 612 d 2 2 - f>22d 2 l) 2 - (f>22C21 - f>12C22) 2 = 0 

(22) 
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cq C '(22) b th rin8 ^ ^ Al ’ ^ A2 arC lhe tW ° ro0tS ° f ° r ’ aftCr Subslitutin 8 e T (27) 


X\ + X 2 = ~ 2 {b\ 2 d 2 \ + b 22 d 22 ) 
and 

■^1^2 = —(612^22 — b 22 d 2 \) 2 
~(b 22 c 2 \ - 612C22) 2 


( 23 ) \\ = - ^ 2 622 d 22 + 2 ^-^22 j 

Ai = — 2 b 22 d 22 ^1 + 

( 24 ) The Cas e of L.E. Control Only 


( 28 ) 


the following conclusions result from eq. (22): 

1 . Because the constant element (independent of X ) 
is always negative (or at best equal to zero), there 
will always be a positive root Ai and a negative 
root A 2 . 

2 . The elements C21 , C22 do not affect the sum (A 1 + 
A 2 ), but they do reduce the constant element, thus 
implying that while they increase the value of Ai , 
they lead to a decrease in A 2 by an equal amount. 

3 . Because for an unbalanced T.E. control, 612 and 

622 are both positive (see appendix A), Ai > -A 2 

if d 2 ] or d 22 or both are negative so as to cause 
(see eq. ( 23 )) 


612^21 + 622^22 < 0 ( 25 ) 

4 . For the absolute value of A 2 to be as small as 
possible while increasing A 1 (to minimize the un- 
stable root and maximize the stable one) without 
changing the control unbalance, one should aim 
at letting (see eq. ( 24 )) 


C21 = C22 = 0 or — = — 

C22 b 22 

with d 2 \ <0, d 22 < 0 ( 26 ) 


5 . 


The optimum control law using the above crite- 
ria for inertially unbalanced T.E. stability is such 
that it satisfies eq. ( 26 ) above and also satisfies 
the equation 


d 22 _ 622 
<^2l 612 

making Aj A 2 = 0 (see eq. ( 24 )) 


( 27 ) 


6. Equations ( 26 ) and ( 27 ) ensure that Ai > 0 (thus 
is stabilizing) while keeping A 2 = 0 (and there- 
fore, A 2 does not contribute to instability). Note 
that A] in this case assumes the value 


^ l — — 2( b\ 2 d 2 \ + b 22 d 22 ) 


Following the above analysis for the case of the 
T.E. control only, a similar analysis will be performed 
for the case of L.E. control only. Hence, in this case, 
one can write 


&12 - 622 = c 2 i = C22 = d 2 \ = d 2 2 = 0 (29) 

Substituting eq. ( 29 ) into eq. (20) the following char- 
acteristic equation is obtained 

A 2 + 2A(6 I1 d 11 + 621^12) - (&ndi2 - b 2 i d u ) 2 

~( b 2 \ C]i — 611 cj2 ) 2 = 0 
( 30 ) 

Equation ( 30 ) is similar in form to eq. (22), and the 

following similar conclusions can be drawn: 

1 . Because the constant element, which is indepen- 
dent of A, is always negative, there will always be 
a positive root A j and a negative root A2 . 

2. The elements cn and C12 do not affect the sum 

+ X 2 , but they do reduce the constant element, 
thus causing an increase in the value of A 1 while 
decreasing A 2 by an equal amount. 

3 . Because an unbalanced L.E. control yields 
(shown in appendix A) positive 61 1 and negative 
621 , the value of Ai will be larger than the abso- 
lute value of A 2 , that is Ai > -A 2 if d \ 2 > 0 or 
•f du <0 (or both conditions), so as to cause 

1 611 dii + b 2 \d\ 2 ] < 0 ( 31 ) 

4 . For the absolute of X 2 to be as small as possi- 
ble, while increasing Aj (to minimize the unsta- 
ble root and maximize the stable one), without 
changing the control unbalance, one should aim 
at letting 


cn = C12 = 0 for — -- — 

\ Cn 6n 

dn < 0, d \ 2 >0 


( 32 ) 


5 . 


The optimum control law based on the above cri- 
teria for inertially unbalanced L.E. control should 
have gains that satisfy eq. ( 32 ) and also satisfy 
the equation 


du _ 6u 
d \2 b 2 \ 


( 33 ) 


6. Equations ( 32 ) and ( 33 ) ensure that \\ > 0 (the 
stabilizing root) while keeping X2 = 0 . There- 
fore, A2 does not contribute to instability. Note 
that \ \ in this case assumes the value 


\\ = — 2 (b\\d\\ + 621^12) 


or, after substituting eq. (33) 


or 


\\ = — 2 b\\d\\ 


2 621 dn -~~ 
0 11 


\\ = — 2 61 1 di 1 |^1 + ( 34 ) 

The Case of L.E.-T.E. Controls 

The characteristic eq. ( 20 ) includes terms associ- 
ated with L.E. control only, terms associated withT.E. 
control only, and also coupling terms between L.E. and 
T.E. controls. If the optimum control law conditions 
derived earlier for the L.E. (alone) system and T.E. 
(alone) system are applied to eq. (20), they ensure that: 

1 . The sum of (\\ + \ 2 ) > 0 since 

(611^11 + b\ 2 d 2 \ + 621 di2 + b 22 d 22 ) < 0 ( 35 ) 


3 . Based on conclusions ( 1 ) and ( 2 ) above, it fol- 
lows that the optimum L.E. (alone) control law, 
and the optimum T.E. (alone) control law, yield 
a combined L.E.-T.E. control law with positive 
Xj and positive \ 2 . Hence, the combined L.E.- 
T.E. control law will always be stabilizing from 
the point of view of inertial unbalance. 

There remains to determine at this stage the com- 
patibility between the unbalanced inertial stability re- 
quirements and those required by the aerodynamic en- 
ergy method to stabilize the system as a result of acting 
aerodynamic forces. 

Summary of Results 

The very interesting results obtained so far are re- 
iterated before proceeding to deal with their compati- 
bility with the aerodynamic control laws: 

1 . Both L.E. alone and T.E. alone mass-unbalanced 
control systems always lead to inertial instabili- 
ties (or neutral stability at best), irrespective of the 
control laws employed. The best control law for 
these single control surface systems may reduce 
these instabilities to neutrally stable oscillations. 

2 . The combined L.E.-T.E. control system is shown 
to be much more powerful than any of its com- 
prising components. With a proper control law, 
the mass-unbalance terms can be made to con- 
tribute to the stability of the system, irrespective 
of its characteristics and its responses. 


2 . All the negative quadratic terms that appear in the 
constant term in eq. (20) vanish, leaving the fol- 
lowing coupling terms: 


COMPATIBILITY BETWEEN 
AERODYNAMIC AND INERTIAL 
CONTROL LAWS 


4 [ 6ll dll £>22^22 + 612^21621^12 ] 

— 2(bndi2 + 621 du) (612^22 + 622 CZ21 ) 

Considering eqs. ( 26 ) and ( 32 ), and remembering 
that 6n, 612, and 622 are all positive, while 621 
is negative, it follows that both terms within the 
square bracket are positive (they add up) and the 
product terms end up being positive, thus adding 
to the square bracketed term. Also, within each 
of the parentheses forming the product, the terms 
add up so that the result is positive. 


Aerodynamic Energy Control Law 

The basic aerodynamic energy L.E.-T.E. control 
law requires the matrix [T] to assume the form 


IT] = 


m = 


<n 1 12 
^21 t22 


an + ien 


022 + 1622 



’ -1 

1 


1 

0.7 


( 36 ) 
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where an terms can assume either positive or negative 
values, and the e„ terms must be positive for aerody- 
namic stabilization (Nissim 1977). The constant terms 
in the matrix in eq. (36) were obtained in Nissim (1977) 
using numerical optimization, maximizing the aerody- 
namic energy eigenvalues over a whole range of sub- 
sonic Mach numbers. Note that when on = e n = 
0 , cq. (36) yields a T.E. system, and similarly, when 
022 = ^22 = 0 , eq. (36) yields an L.E. system. Note 
also that ratios tn /t 12 andf 22 /t 2 i are constant. In the 
following a comparison is made between the aerody- 
namic control law given in eq. (36) and the inertial con- 
trol law given in eq. (17), repeated here for the sake 
of clarity. 


IT] 


cn + tan C12 + idu 
C21 + itfel C22 + id>22 


It is possible to vary these ratios between the t 
terms by changing the first column [—1 lj r in 
cq. (36) relative tr the second column. Such a vari- 
ation may be considered if one is required to sat- 
isfy (or nearly satisfy) the relations in eqs. (24), (25), 
(30), and (31). If this is done, aerodynamic per- 
formance may be somewhat degraded. Therefore, 
relative changes in the first (or second) column of 
eq. (36) can be made during the design stage, when 
the mass unbalance detrimental effects need be over- 
come. In the following, the compatibility between 
eq. (36) and eqs. (24), (25), (30), and (31) is studied 
for the T.E. control system and for the L.E. control 
system. The results pertaining to an L.E.-T.E. con- 
trol system follow the results of the preceeding two 
separate systems. 

Compatibility Between Inertial and 
Aerodynamic Damping for a T.E. System 


Equations (26) and (27) require that 

C22 _ &2_ _ &22_ 

C21 d,2 1 b 12 


(37) 


Eq. (36) indicates that the first equality in eq. (37) is 
always satisfied, that is 


C22 _ ^22 
C 21 <^21 


(38) 


irrespective of the relative values between the two 
columns in eq. (36). As shown in appendix A, 
eq. (A- 8 ) 



(39) 


with both 622 and &12 being positive, for unbalanced 
T.E. control surfaces. The 622 / &12 ratio will be slightly 
larger than 1 for completely unbalanced T.E. control 
surface and can be made to tend to infinity for a stati- 
cally balanced T.E. control (see eq. (A-9)). The stati- 
cally balanced T.E. case is of little interest because in 
this case 612 = 0 , and 622 is very small, thus yielding 
very small numerical values for \ that can be ignored, 
considering the structural damping of the system. 

It follows from the above discussion that the value 
of 0.7 in the second column may need to be increased 
(relative to 1 in the first column) to a value somewhat 
larger than 1. The actual value can be determined in 
the design stage should such a need arise, when con- 
sidering a specific numerical system. It can be readily 
shown that increasing the 0.7 gain moves the single 
sensor for the T.E control system further downstream 
(along the chord) than the 65 percent chord point asso- 
ciated with the 0.7 relative gain (Mukhopadhyay and 
others, 1981). 

In eq. (26) both cfei and dz2 must be negative. 
This requirement corresponds to e 22 assuming nega- 
tive values, and it therefore relates to the case where 
the resulting aerodynamic forces have a detrimental ef- 
fect on stability. Hence, at the region of flutter, aerody- 
namic and inertial dampings are not compatible . This 
is a very important new result. Because flutter occurs 
at relatively low structural frequencies, a design should 
aim at producing a transfer function that yields: 


1 . Large positive values of e 22 around the flutter fre- 
quency (eq. 36) 


2. Negative values for e22 at high frequencies 


3. Value of <222 that tends to zero at high frequen- 
cies (unless eqs. (26) and (27) are completely 
satisfied) 
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4. Value of e22 that tends to zero (with e22 being 3. Value of en that tends to zero (with en being pos- 
negative) at high frequencies, because there is no itive) al high frequencies, because there is no de- 
desire to control the very high frequency modes. sire to control very high frequency modes. 


Compatibility Between Inertial and 
Aerodynamic Damping for an L.E, System 

Equations (32) and (33) require that 

£11 = Jil (40) 

c 12 d\2 021 

Eq. (36) indicates that the first equality in eq. (40) is 
always satisfied, that is 


ci l _ dit 
C12 dl2 


(41) 


irrespective of the relative values between the two 
columns in eq. (36). As shown in appendix A (see for 
example eq. (A- 10)), 


for a mass-unbalanced L.E. control. As the control is 
balanced, the above ratio decreases until it reaches the 
value of zero for mass-balanced L.E. control. Further- 
more, 6n is always positive (or zero when balanced), 
and 621 is always negative, so that 611/621 is always 
negative. From eq. (36) it can be seen that d u /d\ 2 
(and therefore c\\/cn by virtue of eq. (41)) is equal to 
— 1 . If the L.E. control is totally unbalanced, it is pos- 
sible that the 1 in the second column of eq. (36) needs 
to be decreased. If 1 61 1/621 1 < 1 (for a partially bal- 
anced L.E. control), then it is possible the above value 
of 1 needs to be increased. Here again, the actual value 
can be determined later in the design stage, when con- 
sidering a specific numerical system. 

Equation (32) requires that d\ 1 be negative and 
di2 be positive. Equation (36) shows that this require- 
ment is fully met when en is positive. Hence, it can be 
seen that there exists a compatibility between the aero- 
dynamic Ihimphiglmd^ Hence, 

the L.E. control surface transfer function should aim 
to yield: 


1. Positive values of en over a frequency range that 
spans the flutter frequency and extends to very 
high frequencies. 

2. Value of an that tends to zero at high frequen- 
cies (unless eqs. (32) and (33) are completely 
satisfied). 


GENERALIZATION OF RESULTS 
FOR A REALISTIC WING STRIP 


The analysis in appendix A that yields the matrix 
[B c ] essentially assumes a rigid rectangular wing strip, 
performing pitch-plunge oscillations. A more realistic 
representation of the oscillation of a wing strip would 
allow the strip to rotate about its chordwise section, 
while perfonning the plunge motion. In addition, vari- 
ations in both angle of attack and chord lengths along 
the strip always exist and their effects on the [B c ] ma- 
trix need to bc investigated. This latter investigation 
is important because the conclusions reached for con- 
trol laws aimed at coping with mass-unbalanced con- 
trol surfaces depend heavily on both the signs and rel- 
ative sizes of the [ B c ] matrix. 

The analysis of such a generalized wing strip is 
carried out in appendix B. The results yield the fol- 
lowing [B c ] matrix (see eq. (B-12)): 



Ml 

1 + \ r/ W3L + b y yQL 
\' l R/ J 

1 

ZGL ] 

1 

[B e ] = Mb\ 

Ml 

1 + (zifl) WJL + 2 b y ijGL 

X | 


. [~( r G£ + x gl) + ( X L + 0 -4)a:GL 1 2 


M t + \J^ J Vgt + 6j,j/gt 


zgt 


Mt + (5^ VGT + 2 byi/GT 

[( r Ga + s-gt) + ( %t + 0 .4):tgt] 

where M is the mass of the wing strip including the 
control surfaces (see eq. B-10). 


U T =h$ 


(43) 


where Ml and Mt arc the masses of the L.E. and T.E. 
control surfaces, respectively. The parameters h v , a y , 
and by are defined by 


, dh 

n y =-r~= const 

dy 
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* v d(y/b R ) COnSt 

» db 

b v = = const 

dy 

where denotes the strip’s reference semichord 
length. Suffix R relates to the reference section (see 
figure in appendix B), and h R , j/gl. and I/GT denote 

hR = h R /bR 
yoL = yGL/b R 

VGT = IJGT/bR 

The parameters ycL and yGT denote the spanwise dis- 
tances from the reference section of the center of grav- 
ity of the L.E. and T.E. control surfaces, respectively. 

It should be noted that the effect of the generalized 
strip on the [B c ] matrix can be large since h y / J ir and 
oty/ciR can assume large numerical values, with either 
positive or negative signs. It appears, therefore, that 
no single control law can stabilize the generalized strip 
unless the reference chord is chosen such that 

t/GL = i/GT = 0 (44) 

When eq. (44) is satisfied, the expression for [ B c ] 
reduces to the simple 2-D model case except for 
the mass ratio terms. Therefore, all the conclu- 
sions reached for the 2-D pitch-plunge strip ap- 
ply to the generalized strip, provided eq. (44) is 
satisfied. Hence, we reach the following very 
important conclusions: 

1 . The spanwise reference section of an active strip 
should be chosen such that it passes through the 
center of gravity of the control surface. 

2. If bothL.E. and T.E. control surfaces are activated 
along the same strip, these two control surfaces 
should be aligned so that their centers of gravity 
lie along the same spanwise reference section. 

3. If a similar analysis is performed for the aero- 
dynamic forces, based on 2-D aerodynamics, the 
reference section should pass through the center 
of area of the generalized strip. 

4. For best control of both inertial and aerody- 
namic destabilizing forces, the inertial refer- 
ence section described in item (2) above should 
be made to pass through the center of area 


of the wing mentioned in (3) above. Uni- 
form wing and control surfaces do yield such 
a congruency. 

NUMERICAL EXAMPLE 
AND RESULTS 

A mathematical model of the YF-17 aircraft 
(Fig. 1) is used to test some of the results obtained in 
the present work. The available mathematical model 
allows for two activated control surfaces: one leading- 
edge control and one trailing-edge control. (Hwang 
and others, 1978, gives more details relating to the 
mathematical model.) However, no data is available 
regarding the alignment of the center of gravity of L.E. 
and T.E. controls, and therefore, deviations from some 
of the conclusions for the 2-D strip should be expected. 
To keep the focus of this work on the effects of mass 
unbalance on the stability of an active control sys- 
tem, no 11 utter calculations are made, and the system 
is tested for zero dynamic pressure. This means, in 
essence, that eq. (1) is solved for the case where 

[Ar + iAi] = 0 and {F} = 0 

The resulting equation reduces, in effect, to a vibration 
problem involving elastic and inertia forces only, with 
inertial coupling terms between activated control sur- 
faces and structural degrees of freedom. This equation 
can be brought to the canonical form of an eigenvalue 
problem, the solution of which yields the state of sta- 
bility of the system. AH the eigenvalues presented in 
the attached tables relate to the solution of eq. (1) un- 
der the aforementioned conditions. Hence, a negative 
real part of an eigenvalue means a stable structural de- 
gree of freedom, and a positive real part of an eigen- 
value will mean an unstable structural degree of free- 
dom. These eigenvalues should not be confused with 
the energy eigenvalues discussed earlier in this work. 

Table 1 presents the eigenvalues for the 10 elastic 
modes, both for the open-loop case and the closed-loop 
case, using real transfer functions only (that is, using 
real values for t„). It can be seen that a single con- 
trol surface is activated each time, and that the effect 
of this activation on the stability of the system is abso- 
lutely negligible. Table 2 presents similar results per- 
taining to the simultaneous activation of both L.E. and 
T.E. control surfaces. Hence, it can be concluded that 
for the specific example in hand, the real parts of the 
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transfer function have a negligible effect on the stabil- 
ity of the system. Table 3 presents results pertaining to 
the complex transfer function 


o(w„) 


-4 s 2 

s 2 + 4 s + 4 


(s + w„) 


(45) 


where s denotes the Laplace variable. This transfer 
function was chosen to simulate the imaginary compo- 
nent of f,j because the computer program available for 
the solution of eq. (1) is constrained to transfer func- 
tions where the order of the polynomial in s of the nu- 
merator is less or equal to the order of the polynomial 
in s in the denominator. The value of w n was chosen 
to equal 50 to ensure some gain at high frequencies. 
It should be noted that a(u„) introduces the positive 
values of e^. Similarly, — a(w n ) introduces the nega- 
tive values of e,-,-. Contributions of the real part of the 
transfer function are expected to be small, following 
the results presented in Tables 1 and 2. On the basis of 
the foregoing analysis, it is expected that the T.E. with 
— a( 50) and the L.E. with a(50) will each yield the 
most stable system in relative terms, as far as mass- 
unbalance effects are concerned. Table 3 confirms 
the aforementioned expectations. Table 4 presents re- 
sults for a combined L.E.-T.E. active system using the 
a(50) transfer function. Here again, it can be seen 
that the most stable system is the one obtained using 
l n = o(50) and t 2 2 = -o( 50), as predicted herein, 


At this point it should be stressed, once again, that the 
results of the present analysis, based on a general 2-D 
system, were applied to a YF-17 mathematical model 
with no knowledge of the degree of mass unbalance, 
control surfaces’ center of gravity locations, or span- 
wise locations of the control surface centers of gravity. 
Nevertheless, the correlation between analysis and re- 
sults is indeed impressive. 


CONCLUSIONS 

The study of the effects of mass-unbalanced con- 
trol surfaces on the stability of the closed-loop system 
indicates that: 

1. A single L.E. or a single T.E. mass-unbalanced 
control surface may lead to instability (or neu- 
tral stability) arising from inertial forces, irre- 


spective of the control law used, because not 
all energy eigenvalues can be made to assume 
positive values. 

2. The most efficient single L.E. or T.E. mass- 
unbalanced control surface can be made to yield 
a zero energy eigenvalue in addition to a positive 
eigenvalue, thus indicating the possibility of neu- 
tral stability (or an uncontrollable mode). 

3. A combined L.E.-T.E. mass-unbalanced control 
system is shown to permit stabilization, irrespec- 
tive of the mode of oscillation of the wing. 

4. Inertial and aerodynamic control laws can be 
made compatible for an L.E. control system. In- 
compatibility exists, however, between aerody- 
namic and inertial control laws for a T.E. control 
system. 

5. For mass-unbalanced stabilization to be insensi- 
tive to mode of oscillation for a continuous sys- 
tem like a wing surface, the sensor should be 
placed along a streamwise section that passes 
through the spanwise center of gravity of the con- 
trol surface. 

6. When two control surfaces are activated on one 
strip of the wing (like an L.E.-T.E. control sys- 
tem), the two control surfaces should be aligned 
so their spanwise centers of gravity lie on the 
same streamwise section, with the sensors placed 
along this section. 

7. Results obtained herein suggest that the best 
geometrical arrangement for flutter suppression 
should aim at aligning the control surfaces’ span- 
wise centers of gravity (section (6) above) with 
the spanwise centroid of the wing strip, where the 
control surfaces are located. 

8. Numerical results relating to a YF-17 mathemat- 
ical model appear to agree with the theoretical 
predictions based on the analysis made in the 
present work. 

9. Results obtained herein can be used to stabilize 
elastic modes of large space structures by means 
of inertial forces. 
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APPENDIX A 

DERIVATION OF CONTROL SURFACE CROSS-INERTIA TERMS 


s- 30 percent chord point (x = -0.4b) 



9029 

Figure Al. A L.E.-T.E. Control System 


The following derivation of control surface cross inertia terms is performed for a general configuration of a 
L.E.-T.E. system. Numerical estimates, however, will relate to 20 percent chord control surfaces. 

The kinetic energy of the system described in the above sketch is given by: 


K.E. 



(x+ 0 .4 6)a + (xi — x)0] 2 + 
+ (x + 0.4 6)o'+(x — xt) 6] 2 


1 r x T 

— / dm[h + (x + 0.4 6)d] 2 
^ dx L 


(A-l) 


where x L and x T (shown in the figure) denote the x coordinates of the L.E. and T.E. hinge lines, respectively. The 
kinetic energy coupling temis with control surface deflections are given by 

r x L 

(K.E.) C = / dm[(x+ 0.4b)(x L - x)at0 + (x L - x)h0] 

J ~~b 

+ [ dm{(x + 0Ab)(x-x T )a6 + (x - x T ) 'hb] 

J Xt 


or 


(K.E.) C — imizoLhfi + Tfi r pX(y\ hS 


+ 


+ 


f XL 

/ dm(x - x L + x L -f 0 A b)(x i - x)a/3 
J — b 

f dm{ x — xt + xt + 0 A b) ( x — xt) 6ib 

J Xt 


where m L is the mass of the L.E. control surface, m T is the mass of the T.E. control surface, x GL is the distance 
from the hinge line of the center of gravity of the L.E. control, positive if upstream of the hinge line, and x CT is the 
1 1 stance from the hinge line of the center of gravity of the T.E. control, positive if downstream of the hinge line. 

The above equation can be integrated to yield the following form: 


/is i- \ 1 2 ^ n ZGT , i /l i 

•)c - mi - b -0 + mx—^-b —6 + ( Is + (xt + O^^mT-xcTjafS 

+ f — Ip + ( xj, + 0 .4 b) mi xgl ] a/3 


(A-2) 



where 


fO ry 

I& = ( x ~ x t) dm 
J I? 

/ Xl 2 

(x- xt) dm 

and I, and /, are the ntomen. of inertia, about the con, ml surface hinge line of the L.E. and T.E. controls, respective,. 
Equation (A-2) can further be reduced to the following form: 


(K.E.) C = m L XGtb 2 ^P+ m T iGlb 2 -6+ m T b 2 [(rQ 6 + x GT ) + {xt + 0.4)iGT]a<5 


+ m L b 2 [-(rQ H + xgl) + (^t + 0.4 )igl]Q!/5 


(A-3) 


where rot and fo. are the normalized radii of gyration (normalized with respect to the semichord length 6Xof 

^ S l E control surfaces, respectively, about their mspective control surface center of gravity pomB^The 

T indicates that the parameters are normalized with respect to the 
bar above the parameters zgl, ^gt* and x t indicates uiai t 

scmichord length b. 

The expression for the kinetic energy in eq. (A-3) is now differentiated as required by Lagrange’s equations, to 
obtain the inertia terms, that is 

— ( d( K - E -^ - m L XGLb 2 $ + mTXGrb 2 b] 

dt l d(h/b ) J 


± ( a ^m = ,n T i> 

dt \ da J 


(rhs + x gt) + (ir+ 0.4)5 ;gt 


2 

+ rri[ y b 


-(*G/3+ *GL> + (*i + 0.4 )igl 


P 


(A-4) 


Equation (A-4) yields the following coupling mass matrix (with the h/b and a degrees of freedom): 


[Bel = 


mib 2 xci 


mrb 2 igt 


m l b 2 [-{rl j} + x 2 g O + (xl + 0.4)5;gl] m T b 2 [(fl s + x 2 g t ) + (*r+ 0.4)x G t^ 


which can be reduced to 


[ B c ] = mb 2 


mix gl t n o. t 

+ xgl) + (it + 0-4 )igl] rhr[(rG£ + 3 : gt) + + 0.4)x G t 

where m is the mass of the two-dimensional section, m L and m T are given by 

mi 


mrx gt 


(A-5) 


mi = 


mr - 


m 

mT 

TO 


(A-6) 


Note: It should bc noted that if [B c ] is denoted by 


[ B c } 


bn b 12 
£>21 £>22 


(A-7) 
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then cq. (A-5) indicates that 611 and b n are always positive for mass unbalanced L.E. and T.E. control surfaces For 
a 20 percent chord T.E. control, x T = 0.6 thus yielding positive values for b 22 (which will be positive for all TE 

/ C r ° f m T 1CSS ?“? 70 PCrCent chord - Ulat is for > -0 .4 -namely for al] practical cases). Note 
also that & 22 /O 12 will always be larger than 1 , that is 


f>22 

> 1 

012 


( A-8) 


For a 20 percent L.E. control surface x L = -0.6, thus yielding a negative value for b 2X (which will remain 
negative for L.E. control surfaces of sizes less than 30 percent chord-that is for x L < -0 .4 -namely again-for 
all practical cases). J 6 

Numerical Example 

To get some insight into the different values of the 6 t; terms, the following example is presented. Consider 
percent control surfaces (i.e., x L = -0 . 6 , x T = 0 . 6 ) yielding the following values for [B c ], (see eq. (A- 5 )): 


[B c ] = mb 2 ™ L * GL . 2 

-r I /v* “ 


m T XGT 


m L[-(r G/3 + x GL ) - 0.2 zgl] mrt(r Gfi + x^ T ) + z GT ] 


Note that in this case 


^-=\+ th t *GT 

b 12 XGT 


b 2 i —0 .2 XGL — fn n ~ X. 


X GL 

- r 2 r 2 
r Gp~ X GL 


(A- 10) 


Assume also that the airfoil and control surfaces can be approximated to the homogeneous cross section shown below. 

20-percent trailing 
edge control 

I 1 -r^/ 


/ 

20-percent leading 
edge control 

c 

Figure A 2 . Schematic Approximation of L.E.-T.E. System 


Under these conditions, it can be readily shown that 


rn L = 0.22 rh r = 0.11 
XGL = 0.2 XGT =0.133 
X G0 = 0.0133 f^ s = 0.00889 
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Hcncc, 


or 


[B c ] = mb 2 


0.22 x 0.2 

0.22[ -(0.0133 + 0.04) - 0.2 x 0.2] 


0.11 x 0.133 

0.11[0.00889+ 0.0178 + 0.133] 


[ B e ] = mb 2 


0.044 0.01463 

-0.0205 0.01757 



APPENDIX B 

DERIVATION OF CROSS-INERTIA TERMS FOR A GENERALIZED WING 

STRIP WITH CONTROL SURFACES 


Reference 



Figure Bl. Generalized wing strip 

Consider the wing strip shown above. An attempt will now be made to generalize the results obtained in ap- 
pendix A to nonrectangular wing strips which exhibit spanwise deformation in twist, together with rotation (about a 
chordwise axis) which results from the bending of the wing. This implies that the plunge h at the 30 percent chord 
location, the local angle of attack a and the local semichord length b along the strip of the wing can be allowed to 
vary along the span of the strip. It will now be assumed the above variation to be linear and having small derivatives 
(with respect to the spanwise coordinate y). That is, let 

h = Hr + h y y 

a = aR + a y y 

b — bn + b y y ; b y « 1 (B 1) 

where 

dh , . . 

h u = — = const along the stnp 

dy 

a „ = — = const along the strip 
y dy 

fc = — = const along the strip (B-2) 

dy 

and where suffix r relates to the reference section indicated in the sketch above. 

Equation (A-3) can be made to be applicable to an infinitesimal strip. Assume that mi and rar in eq. (A-3) 
now denote local masses per unit span of the L.E. and T.E. control surfaces, respectively. Hence, the infinitesimal 
strip masses will be given by m L dy and m T dy, and eq. (A-3) will now assume the following form (for the spanwise 
integration of the kinetic energy), that is 
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f 2 . . ra 2 

(K.E.) C - / ?^L(y)^GL bh/3dy + / ttit{ y) xqt bhbdy 

* /a i ^01 

rai 

+ Ja mr ^k 2 [(r 2 J(5 + ^gt) + (%t + 0 A)xQj]a8dy 
+ J a rn L (y)b 2 [-( r 2 Gp +x 2 GL ) 

+ ( + 0 .4) xgl ] cv/jdy (B-3) 

Assume also that the local control surfaces chords, relative to the wing chord arc independent of the spanwise 
coordinate y, and that the control surface satisfies the following relations (see appendix A) 


xl = const 
xt = const 
iGL = const 
xgt = const 
rG6 = const 
r G p = const 


(B-4) 


Substituting eq. (B-l) into cq. (B-3) yields 


(K.E.) C 


[J ai m ^y^ b R + b v v)(hR+ h y y)dy XG\j+ |jf 2 m T (y)(b R + b y y)(h R + h y 

* [f 


y)dy 


xqt 6 


a\ 

&2 


m T (y)(b R + b v y) 2 (a R + a y y)dy J [(r ( 2 56 + z£ T ) + (i r + 0.4 )x G t] 6 

J ai m ^y^ b R + b„y) 2 (a It + a y y)d^ \-(f 2 Gf) + xh ,) + (x L + 0 .4)i GL j /? 


XGI.P 


Expanding eq. (B-5) and ignoring second order terms in b y , one obtains 

/ ° 2 

i m L(y)(b R h R + b R h y y + b y h R y + b v h y y 2 )dy 

/ a>2 

m T (y)(b R h R + b R h y y + b y h R y 4 b y h y y 2 ) dy\x G1 8 

r f ai 

+ / rn T (y)(b 2 R + 2b R b y y) (a R + ay) dy 

.J a i 


x t ( r G6 + ^gt) + (xt + 0 A)xqx]8 

fQ2 

+ / + 2 b R b v y){a R + d v j/)dj/ 

x ( — ( *G0 + ^gl) + (xl + 0 .4)i GI ]/? 


(B-5) 


which can be reduced to 

(K.E.) C = [M L b R h R + M L b R h y y GL + M L b y y CL h R + M L r 2 yL b y h y }x GL p 

+ [ M T b R h R + M T b R h y yc, T + M T b y y Gy h R + M T r 2 yT b y h y ]x Q1 6 
+ [ M T b 2 R a R + M T b 2 R a y y GT + 2 M T b R b y a R y GT + 2 M T b R b y a y r 2 yT } 

* [(r Gfi + *ar) + ( X T + 0 .4).x cl ] 6 + [M L b 2 R a R + M L b 2 R a y y GL + 2 M L b R b y a R y Gh + 2 M L b R b y a y rl L \ 

x r-(r G/3 + 5 G l) + (ii + 0.4)i GI . ]/3 , R , 
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where Mr and Mr denote the masses of the L.E. and T.E. control surfaces, respectively, and r, L .r, T *“ttlhc 
Tadh onnenia about die reference streamwise section of die L.E, and T.E. conuol surfaces, respecuvely. Vfldiou 
affecting the generality of the results, assume that a relationship ex.sts between h y and h R , a y and a R ( 
relationship existing during a modal oscillation), such that 


hy_ _ 

h r hft 

Oy_ = Oy_ 

oiR a R 


so that cq. (B-6) can be written in the form 
( K.E.) C = M L b 2 R 

+ Mt^r 


1 + (fe) ^ + (fe) ^ + ( ^ ) (£)] 

*) * OT * (fe)] iGT (j 


1 + 


h R J yGT + \b R , 


' XGL (^) ^ 


h R 


+ Mrf.fi |^1 + (^) HOT + 2 (^) KT + 2b s b, (f£) ) j 

X [(fofi + igt) + (zr + 0 - 4 )xg t ] a R S 
+ Ml&[l + (f*) ML * 2 (^) yon + !*A (^) (-jf) 

x [— (f'G/? + ^gl) + (%l + °- 4 )^gl] Oi R $ 


Assuming 




(B-7) 


(B-8) 


\ w Ilr r 

one can neglect the second order terms in eq. (B-8) involving b v ) and b y ) • Differentiating eq. (B 8), as 
specified by Lagrange’s equations (see appendix A), one obtains the following [ B c 1 matnx 

. _ < . v / i \ 


[B c ] = Mb 2 R 


Ml [l + (fe) l/GL + (fe) vgl] XGL 


L \"«/ \-xw -» I 2 

al [l + (a) VOL + 2 (jj) yc L ] X j Mr[l + fe) Vcr + 2 (&) 

1 I r . <u / - r\ A\ : 


L ' “ ' ' 

f-(rGfl +i GL) + (^ + 0 - 4 > iGL ] 


*r[l + (fe) vex + (fe) ^gt] Sgt 


I T/GT X 


Uric + x 2 m) + (*r + 0 . 4 )xgt] 


(B-9) 


where M is the mass of the strip and 


M l = 
Mt = 


Ml 

M 

Mt 

~M 


(B-10) 


The element inside the square matrix can be nondimensionalized by normalizing with respect to the reference chord, 
in the following manner: 
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Let 


da 

* V ~dGJbR) 

h.R = /ifi/bs 
1/GL = |/GL /t>R 
j/GT = J/GT /f>R 


Then eq. (B-9) can be written in the following form: 


[B c ] = Mb\ 


Ml 

1+ 1 

(ft) + Wol ] 

Ml 

F 

IS) ito i- + 2 M GL ] 


-<Tg/3 + ®GL) + ( X L + 0- 4 )iCL 


| Mt 

\ 

1 + yGT + M GT 

xgt 

j Mt 

l 1 + (S) ^ GT + 2 M GT 

|x 

| [(r^ + ^gt* + + 0.4)xgt] 


(B- 11) 


(B-12) 
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TABLE 1. Variation With Single Control Mass Unbalance of 
Complex Eigenvalues \ Using Real Transfer Function 


X Real 

R10 

^Imaginary 

Rll 

*ReaI 

^Imaginary 

R12 

*Rcal 

^“Imaginary 

R13 

*Real 

^Imaginary 

R14 

A Real 

^Imaginary 

— 0.1563E— 12 

0.3893E + 03 

— 0.4707E— 12 

0.3884E+03 

0.1066E— 12 

0.3906E+03 

— 0.2132E— 13 

0.389 9E+ 03 

— 0.3055E— 12 

0.3888E+03 

-0. 1421 E- 12 

0.3168E + 03 

— 0.2842 E— 1 3 

0.3094E+ 03 

0.1847E— 12 

0J 249E+03 

— 0.5365E— 12 

0.321 3E+ 03 

0.7105E— 13 

0.3129E+03 

-0.U37E- 12 

0.267 1E+ 03 

— 0.3197E— 13 

02668E+03 

0.4263E— 13 

0 2673E+03 

0.3104E— 12 

02691E+03 

-0.9948 E- 13 

02649E+03 

— 0.3 126 H— 12 

0.2277E+03 

— 0.1616E— 12 

02277E+03 

0.1172E— 12 

02277 E+ 03 

0.2034E— 12 

02278E+ 03 

025 58 E— 12 

02276F.+03 

— 0.2132E— 12 

0.161 3E+ 03 

— 0.1279E— 12 

0.1613E+03 

0.525 8E- 12 

0.1613E+03 

0.2842E— 13 

0.1613E+ 03 

— 0.56 84 E— 13 

0.161 3E+03 

0.1705E— 12 

0.1203E+03 

0.1030E— 12 

0.1203E+03 

02416E-12 

0.1203E+03 

— 0.4832E— 12 

0.1203E+ 03 

— 0.1634E— 12 

0.1203E+03 

-0.341 IE- 12 

0.9351E+02 

0.1151E— 11 

0.9319E+02 

03439E-11 

0.9384E+02 

0.3979E— 12 

0.9357E+02 

— 0.6928E— 12 

0.9345E+02 

0.9948E— 13 

0.8416E+02 

0.126 5E— 11 

0.8404E+02 

-0.547 IE- 12 

0.8427E+02 

— 0.3268E— 11 

0.841 9E+ 02 

0272 IE- 11 

0.8413E+02 

0.1705 E— 12 

0.451 5E+ 02 

0.4062E— 11 

0.4514E+02 

0.5826E— 12 

0.451 6E+ 02 

— 0.4952E— 11 

0.451 7E+ 02 

0.6297E— 12 

0.4513E+02 

0.2842E— 13 

0.2908E+02 

—0.591 2E— 11 

02903E+02 

— 0.4743E— 1 1 

02913E+02 

— 0.8630E— 1 1 

0.2909E+ 02 

— 0.9603E— 11 

02907E+02 


RIO— Open I>oop 
R11-T.E. Only; = 4 

R12-T.E. Only; — —4 

R13-LJL Only;tjj =4 
R1 4-L.E. Only; t j j — — 4 
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TABLE 2. Variation With L.E.-T.E. Mass Unbalance of 
Complex Eigenvalues A Using Real Transfer Function 


A RcaJ 

^Imaginary 

A Real 

^Imaginary 

*Rca 1 

^Imaginary 

^Rcal 

V. 

0.1492E— 12 

0.3889E+03 

-0.2700E-12 

0.3880E+ 03 

— 0.2842E— 13 

0.391 4E+03 

03542E— 1 2 

03899E+03 

— 0.46 90 E— 12 

0.3140E+03 

— 0.4263E— 13 

0305 5E+ 03 

0.7816E— 12 

0.3292E+03 

— 0.5087E— 11 

0.321 IE+ 03 

— 0.4050E— 12 

0.2687E+03 

— 0.1137E— 12 

02648E+ 03 

0.2842E— 13 

0.2695E+03 

— 0.6928E— 12 

03650E+ 03 

0.3979 E— 12 

0.2278E+03 

0.6182E— 12 

03277E+03 

— 0.7105E— 13 

03277E+03 

-02380E— 12 

02276E+03 

-02387E- 1 1 

0.1613E+03 

— 0.8100E— 12 

0.1613E+ 03 

0.7319E— 12 

0.161 3E+03 

— 0.3979E— 1 2 

0.16I3E+03 

— 0.1279 E— 12 

0.1203E+03 

03527E-13 

0.1203E+03 

0.1243E— 11 

0.120 3E+03 

— 03116E— 12 

0.1203E+03 

— 02984E— 12 

0.9325E+02 

-02540E- 12 

0.931 3E+ 02 

— 0.4849E— 12 

0.939 IE+02 

03126E— 12 

0.9377E+ 02 

0.1580E— 10 

0.8407E+02 

0.5329E— 1 1 

0.8401E+02 

— 0.1023E— 10 

0.8430E+02 

03323 E— 10 

0.8425E+02 

— 0.1696E— 11 

0.4517E+02 

0.4245E— 1 1 

0.4512E+ 02 

— 0.3924E— 12 

0.451 HE+ 02 

0.5677E— 11 

0.451 3E+ 02 

— 0.27 56 E— 10 

0.2904E+02 

0.3519E— 11 

039O2E+02 

— 0.1802E— 10 

0 7914E+02 

-0.621 0E- 11 

03912E+02 


RI5-L£.-TJB.;t n = 4, in =4 
R16-LJE.-T.E.; tj j = -4 *22 “4 
R17-L£.-TJE.;( n =4,^ = -4 
R18-L.E.-T B.;f n =-4^t 22 = -_ 4 


TABLE 3. Variation With Single Control Mass Unbalance of 
Complex Eigenvalues A Using Complex Transfer Function 


A Rral 




gtnary 


0.1333E+00 
0.1190E+01 
0.3996E— 01 
— 0.9477E— 02 
0.278 8E- 02 
— 0.1761E— 03 
0.I332E+00 
0.4911 E-01 
0.227 6E- 02 
0.1686E— 01 

R50-TE, Only, t ')') 
R51— T.E. Only; tjj 
R52 — Lit, Only; t-i i 
R53 — LE. Only; tjj 
o(w) is defined in eq. 


RSI 

^Real ^Inru 


lagmary 


0.3893E+ 03 

0.3170E+ 03 

0.267 1E+ 03 

0.2277E+ 03 

0.1613E+03 

0.1203E+03 

0.9359E+ 02 

0.841 9E+02 

0.4515E+02 

0.2912E+02 

= o(50) 

= — o( 50) 

= o(50) 

* — o( 50) 

(45) 


— 0.1318E+00 
—0.1 188E+01 
-0.406 2E- 01 
0.9 56 0E— 02 

— 0.278 9E- 02 
0.1741E— 03 

— 0.1313E+00 

— 0.5072E— 01 

— 02266E— 02 
-0.1678E— 01 


R52 

^ Imaginary 


R53 

A Real X 


0.3893E+03 
0.3166E+ 03 
0.267 1E+ 03 
02277E+03 
0.1613E+03 
0.1203E+Q3 
0.9343E+ 02 
0.841 3E+ 02 
0.451 5E+ 02 
02904E+Q2 


— 0.6984E— 01 
— 0.6363E+00 
— 0.3842E+00 
— 0.1189E— 01 
— 0.8265E— 02 
0. 1721 E— 03 
— 0.2617E— 01 
— 0.1138E— 01 
— 0.9934E— 02 
— 0.2896E— 02 


0.3893E+ 03 
0.3167E+03 
02.67 0E+ 03 
02277E+03 
0.161 3E+ 03 
0.1 203E+03 
0.9349E+02 
0.841 5E+ 02 
0.451 4E+ 02 
02907E+02 


0.7040E— 01 
0.6400E+00 
03808E+ 00 
0.1173E— 01 
0.8219E— 02 
— 0.1718E— 03 
0.261 9E- 01 
0.1129E— 01 
0.9931E— 02 
0.2888E— 02 


0J893E+O3 

0.3169E+O3 

02672E+03 

0.2277E+03 

0.1613E+ 03 

0.1203E+03 

0.9353E+Q2 

0.8417E+02 

0.4516E+02 

0.2909E+02 



TABLE 4. Variation With L.E.-T.E. Mass Unbalance of 
Complex Eigenvalues \ Using Complex Transfer Function 



^Real 

— 0.2010E+00 

^Imaginary 

0.3893E+03 

— 0.1828E+01 

0.3165E+03 

-0.4229E+00 

0.26706+03 

— 0.2363E— 02 

02.277E+ 03 

— 0.1102E— 01 

0.1613E+03 

0.3482E— 03 

0.1203E+03 

— 0.1568E+00 

0.9342E+02 

— 0.6237E— 01 

0.8412E+02 

— 0.1222E— 01 
— 0.1964E— 01 

0.451 3E+ 02 
0.2903E+02 


>Rcal ^Imaginary 


02042E+00 

0J893E+03 

0.1827E+Q1 

0J171E+ 03 

0.4226E+ 00 

02.672E+03 

0.221 3E— 02 

02277E+03 

0.1104E— 01 

0.161 3E+ 03 

— 0.34 60 E— 03 

0.1203E+03 

0.1601E+00 

0.9360E+Q2 

0.6014E— 01 

0.8420E+02 

0.1219E— 01 

0.4517E+02 

0.1978E— 01 

0.2913E+02 


>>Rcal * Imaginary 


0.6289 E— 01 

0.3893E+03 

0.5563E+00 

0.3169E+03 

— 0.3462E+00 

0.267 0E+03 

-02133E— 01 

0.2277E+03 

— 0.551 4E— 02 

0.1613E+03 

— 0.5942E— 05 

0.1203E+03 

0.1064E+00 

0.9357E+02 

0.3798E— 01 

0.8419E+02 

— 0.7639 E— 02 

0.451 4E+ 02 

0.1393E— 01 

02911E+02 


^Real 

Imaginary 

— 0.6206 E— 01 

0.3893E+03 

-03454E+00 

0.3167E+03 

0.3385E+ 00 

0J672E+03 

02133E-01 

02277E+03 

0.5393E— 02 

0.1613E+03 

0.3435E— 06 

0.1203E+03 

— 0.1058E+00 

0.9345E+02 

— 0.3918E— 01 

0.8413E+02 

0.7684E— 02 

0.4516E+02 

— 0.1 393 E— 01 

0.2905E+02 


R54-L£.-T.E.;* n =0(50),^ = -<*(5°) 
R55-L£.-TE.;*]l = -a(50) ,(22 - a(S0) 
R56-L£.-T£.;t u = a(50 ),t n **(50) 

R57-L£.-T£.;tn = -»(50),t22 =-“(50) 
o(w) is defined ineq. (45) 
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